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Abstract 

Context. Large scale motions in galaxies (supernovae explosions, galaxy collisions, galactic shear ...) generate turbulence, 
allowing a fraction of the available kinetic energy to cascade down to small scales before being dissipated. 
Aims. We establish and quantify the diagnostics of turbulence dissipation in mildly irradiated diffuse gas in the specific 
context of shock structures. 

Methods. We incorporate the basic physics of photon-dominated regions into a state-of-the-art steady-state shock code. 
We examine the chemical and emission properties of mildly irradiated (Go=l) magnetised shocks in diffuse media 
(riH = 10 2 to 10 4 cm -3 ) at low to moderate velocities (from 3 to 40 km.s -1 ). 

Results. The formation of some molecules relies on endoergic reactions. In J-shocks, their abundances are enhanced by 
several orders of magnitude for shock velocities as low as 7 km.s -1 . Otherwise most chemical properties of J-type shocks 
vary over less than an order of magnitude between velocities from about 7 to about 30 km.s -1 , where H2 dissociation 
sets in. C-type shocks display a more gradual molecular enhancement as the shock velocity increases. 
We quantify the energy flux budget (fluxes of kinetic, radiated and magnetic energies) with emphasis on the main 
cooling lines of the cold interstellar medium. Their sensitivity to shock velocity is such that it allows observations to 
constrain statistical distributions of shock velocities. 

We fit various probability distribution functions (PDFs) of shock velocities to spectroscopic observations of the galaxy- 
wide shock in Stephan's Quintet (SQ) and of a Galactic line of sight sampling diffuse molecular gas in Chamaeleon. In 
both cases, low velocities bear the greatest statistical weight and the PDF is consistent with a bimodal distribution. In 
the very low velocity shocks (below 5 km.s -1 ), dissipation is due to ion-neutral friction which powers H2 low energy 
transitions and atomic lines. In moderate velocity shocks (20 km.s -1 and above), the dissipation is due to viscous 
heating and accounts for most of the molecular emission. In our interpretation a significant fraction of the gas on the 
line of sight is shocked (from 4% to 66%). For example, C + emission may trace shocks in UV irradiated gas where C + 
is the dominant carbon species. 

Conclusions. Low and moderate velocity shocks are important in shaping the chemical composition and excitation state 
of the interstellar gas. This allows to probe the statistical distribution of shock velocities in interstellar turbulence. 

Key words, shock waves - astrochemistry - ISM: abundances - ISM: kinematics and dynamics - ISM: molecules 

1. Introduction citation diagram in SQ implies a distribution of 

tempatures much above the equilibrium tempera- 
Bulk motions in galaxies generated by e.g. su- ture set by UV and cosmic ray heating. The same 
H ■ pernovae explosions, galaxy collisions and galac- holds on much smaller scales in the case of the dif- 
03 tic shear drive turbulence in the cold interstellar fuse I SM in the Solar N eighbourhood (|Grv et al.l 
medium (ISM). A fraction of the available kinetic 120021 : llngalls et aU 120111 ). The range of gas tem- 
energy cascades down to smaller scales and lower peratures can only be accounted for if the dis- 
velocities. This is spectacularly illustrated by the sipatio n heats the gas thro ugh spatially localised 
observations of the galaxy collision in Stephan's events (jFalgarone et al.l l2005h. This idea is also sup- 
Quintet (SQ). The relative velocities of the galaxies ported by observations of chemical species in the 
(on the order of 1000 km.s" 1 ) would be expected diffuse ISM such as CH+ and SH + which cannot 
to dissipate in high velocity shocks, thus creating a be reproduced by UV alone (jNehme et al.l 120081 : 
warm and hot plasma devoid of molecules. However, iGodard et al.ll2012f). 

one observes th at H2 cooling is greater than the X- These dissipation processes have yet to be stud- 
ray l uminosity (Apple ton et aLl 120061 ; ICluver et alJ ied. Dissipative structures could take the form of 
2010). This demonstrates that energy dissipation e.g. shocks, vortices, current sheets, shear layers 
in the interstellar space involves an energy cascade and most likely involve the magnetic field. This 
and molecular gas cooling. Moreover, the H2 ex- broad variety of processes, their different time 
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scales and therefore different impacts onto the 
global gas energetics and chemistry, call for simpli- 
fied models to quantify observational diagnostics. 
Spectroscopy of molecular H2 with Copernicus in 
the diffuse ISM motivated the first MHD shock 
models propagating in diffuse gas. These first mod- 
els were also key to prov ide a first interpretation of 
CH+ in the diffuse ISM (lMullanlll97ll: iDraine et al.l 
1983|; iFlower et ail Il986l; iGredel et aU 120021) . ISO- 
SWS observations (Falgar one et aLI [^005) of the 
pure rotational lines of H2 in diffuse gas were also 
interpreted in the framework of mild turbulent dis- 
sipation (low velocity MHD shocks and/or veloc- 
ity shears in s mall scale vortices). More recently, 
iGodard et~al . (2009) have modeled turbulent dissi- 
pation bursts in the diffuse ISM as vortices at very 
small scales and computed the chemical signatures 
they imprint on the gas. 

In the present work, we use an updated 
version of the m agnetised sho ck models of 
IFlower &; Pineau des Foretsl (|2003h to quantify 
composition and cooling lines for a range of shock 
velocities (from 3 to 40 km.s x ) including very low 
velocities which were not thoroughly considered 
before. We compute grids of such shocks for two 
strengths of magnetic field and for three different 
densities. Our preshock conditions are representa- 
tive of the cold diffuse ISM. Hence, we have in- 
corporated the treatment of mild irradiation which 
includes UV heating but also impacts the chemistry 
through photo-ionisation and photo-dissociation. 

The shock velocity can be a crucial parameter 
even at very low velocity. As it happens, the gener- 
ation of molecules in the interstellar medium (ISM) 
is initiated by the formation of the H2 molecule on 
dust grains surfaces. Then, the formation of more 
elaborate molecules relies on two main paths. First, 
one can add a proton to H2 to form Hjj~ and then 
transfer the extra proton to a single atom (C, O, S 
and Si, with the notable exception of N). But H^ 
is usually much less abundant than H2. Second, one 
can directly exchange an atom with one proton of 
the H2 molecule. But these reactions are subject 
to energy barriers with high characteristic temper- 
atures (such as 2980K in the case of the oxygen 
atom). Hence both paths are difficult and the ob- 
served molecular complexity of the ISM stands as 
a puzzle. Nevertheless, temperatures required to 
open the second path are obtained in low veloc- 
ity shocks (7.5 km.s -1 is enough to reach 3000K). 
These shocks may therefore play an important role 
in shaping the molecular chemistry of the ISM. 

The results of the present models may contribute 
to the interpretation of the observations of the full 
set of cooling lines of the cold neutral medium (CII, 
CI, OI, H2, CO and H2O) in galaxies observed with 
the Spitzer and Herschel space telescopes. To illus- 
trate this, we build models from statistical distri- 
butions of shock velocities which we compare to ob- 
served data i n the SQ and in a line of sight in the 
Chamaeleon (|Grv et al.l 1200 2 h We use our best fit 
models to discuss the impact of kinetic energy dis- 
sipation on the physical state and chemistry of the 
molecular gas in these two sources. We also make 



predictions on other diagnostic lines to be observed 
with Herschel. We present the numerical method for our 
shock models in section 2, with the properties of our grid 
of shocks in section 3. In section 4, we use our results to 
interpret the H2-excitation diagram observed in SQ and 
Chamaeleon. We summarize and discuss our results in sec- 
tion 5. 

2. Numerical method 

The models we present in this paper are based on the 
plane-parallel steady-state shock code implemented in 
IFlower k Pineau des Foretsl ([2003). We work from their 
version and include few more refinements mainly in order 
to deal with moderate irradiation. 

2.1. Radiation 

We start from the re ac tions network used in 
IFlower k Pineau des Foretsl (|2003l ) and include the 
relevant photo-reactions. The photo-reactions have their 
rates of the form 

R = ctG e- pAv (1) 

where a and ft are constants. We use an incident fiel d equal 
to th e standard interstellar radiation field (ISRF, iDraind 
Il978f ). thus taking Go = 1. The extinction A v is integrated 
along the model from the pre-shock where we use a value 
of A v — A v q =0.1 which models extinction from the irra- 
dation source by a "buffer" of matter. The value of A v o is a 
parameter of the problem which selects how much mass is 
contained in this buffer. The local value for the extinction 
is then computed as 



A 1} 



dz <7 B nH 



(2) 



where er g = 5.34 x 10 -22 cm 3 pc -1 is the effective ex- 
tinction per H nuclei column-density, nn is the lo- 
cal density of H nuclei, x is the current position and 
Xo is the position at the entrance of the preshock (ie: the 
point where we start our simulations). Our shocks are hence 
assumed to be irradiated "backward" compared to their di- 
rection of propagation. However this matters only a lit- 
tle because the total extinction through these shocks is 
small (on the order of AA V = 0.01 (tih/100 cm -3 ), where 
nn is the pre-shock density). This setup is in fact ex - 
actly similar to the one used by Ber gin et al.l (|2004). 
However they focus on the properties of a uniform 
pressure photon-dominated region (PDR) following 
the shock, whereas we stop our computation right 
behind the shock, and they consider pre-shock den- 
sities of nn = 1 cm 3 much lower than ours. 

The reaction rates for the photo-dissociation of H2 and 
CO include an additional factor /shield to accou nt for the 
shield ing and self-s hielding. We use the tabl es of iLee et al.l 
(fl996l ) for CO and lDraine ■ k Bertoldil (fl996l ) for H 2 with a 
Doppler parameter 



b D = 



''turb 



(3) 



where u t h is the thermal velocity of H2 molecules and 
Wturb = 1 km.s -1 accounts for microturbulence. However, 
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we neglect the effect of velocity gradients on the self- 
shielding. For instance: 



Parameter 



Value 



Comment 



where 

/shield 



(2.54 x KT 11 *- 1 ) f m 



-6.3 A v 



0.965 



0.035 e~ 



..5 XlCT 



(i + x/b D y 



\/l + x 



(4) 



(5) 



with x = iV(H 2 )/(5 x 10 14 cnT 2 ). Here, the quantity iV(H 2 ) 
is the column-density of H 2 molecules, computed as for the 
extinction: 



A(H 2 ) = AT (H 2 ) 



dzn(R 2 ) 



(6) 



where n(H 2 ) is the local density of H 2 molecules and iVo(H 2 ) 
is the quantity of H 2 molecules in the "buffer" which shields 
the medium from the external radiation field. We match 
A (H 2 ) ~ 10 20 cm -2 to A v0 = 0.1 by assuming this buffer 
is mainly molecular. Similarly we define ATq(CO) for the CO 
self-shielding but we use A*o(CO)=0 as A v q — 0.1 is usually 
well below the extinction at which CO is self-shielded in 
standard models of photon-dominated regions (PDR) . 

With this irradiation model, we integrated the chemi- 
cal and thermal evolution of a fluid parcel with fixed den- 
sity moving away from the radiation source at a constant 
velocity updr where the subscript PDR refers to photon- 
dominated region. The resulting spatial profile (the trace 
of the thermal and chemical composition history of this 
fluid parcel along its path) corresponds to a steady PDR 
front moving towards the irradiation source at speed fpDR- 
For very small velocities we recover a steady PDR front 
which compares quite sat isfactorily to the Meudon PDR 
code (|Le Petit et al.l I2006T ) for G = 1. In the frame of 
this comparison we switched off grains and PAHs reac- 
tions in order to remain closer to the Meudon PDR 
code setup. We note the PD R structure is insensi t ive to 
the choice of iLee et al.l (|l996h or iDraine" fc Bertoldil (|1996h 
for H 2 photo-dissociation. This is due to the fact that the 
rates differ only in regions where the absolute value of the 
rates are small, hence unimportant. For the same reason, 
we expect our results should poorly depend on the ex- 
act value of the above ment ioned Doppler parameter bo- 
iBlack fc van Dishoeckl (|l987f ) also shows that this parame- 
ter is not crucial in PDRs. 

Compared to IFlower" fc Pineau des Foretsl (|2003l ) , we re- 
freshed the co llision rates of 01 by H ato ms with the rates 
computed by lAbrahamsson et al.l (|2007f h These rates en- 
ter the computation of the atomic cooling due to O atoms. 
Atomic and mole cular cooling are otherwise i dentic al to 
those employed in lFlowe~ fc Pineau des Foretsl (|2003ft . For 
instance, we use the mole cular cooling rates tabulated in 
iNeufeld fc Kaufman! (11993T) for the cooling by CO and H 2 0. 

2.2. ^-excitation 

We follow the time-d ependent excitation of H 2 al o ng th e 
shock structure as in IFlower fe Pineau des Foretsl (|2003h . 
Here, we include the treatment of the H 2 population for 
the 149 lowest energy levels. We check that this allows to 
compute the H 2 cooling accurately for shock velocities up 
to at least 40 km.s -1 for the range of densities we consi der 
in this work (see also IFlower fc Pineau" des Foretsll200l . 



nu 
u 

b = B/n% 



1/2 



A v0 
iV (H 2 ) 
iV (CO) 

Go 

c 

OPR 



10 2 , 10 a or W'cm-' 6 
3 to 40 km.s" 1 
0.1 or 1 

0.1 
10 20 cm- 2 

1 



H nuclei pre-shock density 
Shock velocity 
Dimensionless magnetic field 
Buffer extinction 
Buffer H2 column-density 
Buffer CO column-density 
External radiation field 
Nominal cosmic rays flux 
Preshock H2 ortho/para ratio 



Table 1. Values for the main physical parameters in our 
models. 



We also use the knowledge of the population of the 8 
lowest H 2 rotational levels to compute more accurately the 
rate of the reac tion C + + H 2 . To c omput e the rate state 
by st ate, we use iGerlich et al.l (|l987l ) (as in lAgundez et al.l 
120101 line 2 of ta ble 1) for H 2 leve ls with J = P.. .7 and 
0, and we use iHierl et al.l (fl997h (as in Agundez et al. 



(|2010h . line 3 of table 1) for the other levels 



2.3. Grains 



As inE ower fc Pineau des Foretsl (|2003f ). we treat adsorp- 
tion, collisional sputtering and collisional desorption of 
molecules from and onto grains. However, unlike them, we 
take pre -shock conditions without ice mant les. 

As inH ower fc Pineau des Foretsl (|2003l ) we account for 
the charge of grains by including all charge exchanges in- 
volving grains and electron detachment by cosmic ray in- 
duced secondary photons. The heating due to the photo- 
electric effect is included, however we discard the detach- 
ment of electrons due to the radiation field in the chemical 
network. 



3. Grid of Models 

3.1. Numerical protocole 

For each value of our parameters described in the next sec- 
tion, we integr ate the steady-state equations of mult i-fluid 
MHD shocks (jFlower et all 119851: iHeck et al.l fl990h from 
entrance conditions at thermal and chemical equilibrium 
(see details of the pre-shock chemical conditions below). 
As stated above, we account for the changes in the proper- 
ties of the irradiation field due to the increasing absorption 
as we penetrate deeper into the shock structure. The com- 
putation is stopped when the temperature decreases back 
to 20% above the temperature of the pre-shock. 

3.2. Choice of parameters 

We tune most of our parameters to the typical conditions 
encountered in the dilute interstellar gas in our galaxy. We 
list the main physical parameters of our model in table 13.21 
We specify indirectly the strength of the magnetic 
field transverse to the shock speed by using the non- 
dimensional param eter b — (B/lfiG)/ ^/rtn/cm 3 which 
ICrutcher et al.l (|2010f ) observe to take values from b = 0.1 
to b = 1 for our range of densities. We later refer to these 
values by highly magnetized shocks for 6=1 and weakly 
magnetized shocks for b = 0.1. 
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Fractional value n(X)/nu 



Element X 


total 


gas 


grain cores 


PAHs 


u 
n 


i 
i 


i 
i 






He 


n 1 

U. 1 


U. 1 






c 


3.56(-4) 


1.38(-4) 


1.63(-4) 


5.40(-5) 





4.42 (-4) 


3.02 (-4) 


1.40 (-4) 




N 


7.94(-5) 


7.94(-5) 






Mg 


3.70(-5) 




3.70(-5) 




Si 


3.70(-5) 


3.37(-6) 


3.37(-5) 




S 


1.86(-5) 


1.86(-5) 






Fe 


3.23(-5) 


1.50(-8) 


3.23(-5) 





Tabl e 2. Elemental co mposition of the pre-shock gas (as 
in iGodard et alJ [2009) . Initial conditions assume no ice 
mantles (bare cores). Note the gas to dust mass ratio is 
/°gas/P(fust — 1^0- The numbers in the parentheses denote 
the powers of ten. 



For both of these assumptions on the magnetic field, 
we build three grids of models, one for each pre-shock 
density between nn = 10 2 cm -3 , jih = 10 3 cm -3 and 
7iH = 10 4 cm" 3 . In the following, we will use the term 
low velocity for shocks below 20 km.s -1 , moderate 
velocity above 20 km.s 1 and high velocity at and 
above 40 km.s" 1 . Each grid spans a range of velocities 
from 3 up to 40 km.s" 1 with a step of 0.5 km.s -1 , for an 
integration time of about 6 hours per grid on a typical work- 
station. We choose our minimum velocity of 3 km.s -1 above 
the Alfven speed in the neutral gas in the 6=1 case: 



B 



VA 



Att(^ h / a.m.n.) 



= 1.85 6 km/s (7) 



where the mean mass per H nucleus is p,jj = p/nn = 1.4 
a.m.u. . We choose the upper bound of 40 km.s -1 because 
of the limitations of our shock mo delfO . 

At each fixed pre-shock density, we first evolve the gas 
chemically and thermally during 10 7 years, which brings the 
gas close to thermal and chemical equilibri um. We start 
with t he same elemental abundances as in IGodard et al.1 
(2009), which are summarised in table [U Note the gas 
phase abundance in the pre-shock is set to a tenth of the Si 
l ocked in the gra ins cores, consistent with observations 
(jjenkinsl l2009h . We also include PAHs with a fraction 
n(PAH)/nH = 10 -6 of the representative species C54H18. 
Without irradiation (Go = 0), PAHs influence the ionisa- 
tion degree of the gas and the charge-ion coupling. However 
with mild irradiation (at Go = 1) it turns out that C + ions 
dominate the charge fluid and the role of PAHs is neg- 
ligible. With our current irradiation parameters (Go, ex- 
tinction and H2 buffer), the atomic hydrogen fraction is 
n(H)/n H = 7.9(-2), 1.3(-2) and 2.0(-3) in the preshock gas 
with respective densities njj = 10 2 cm -3 , 10 3 cm -3 and 
10 4 cm -3 . We then use the resulting state as our pre-shock 
conditions to run the shock model for each velocity in turn. 

3.3. C- and J-type shocks 

Steady-state magnetized shocks in the interstellar medium 
are of one of two kinds: J-type shocks in which the ki- 



Indeed, the code does not treat doubly ionised species 
and radiative emission from the post-shock itself which 
might affect pre-shock conditi ons at high densities (see 
iFlower fe Pineau des Foretsll2010t ). 



netic energy is dissipated viscously in a very sharp velocity 
jump (hence 'J') followed by a thermal and chemical re- 
laxation layer and C-type shocks in which kinetic energy is 
continuously (hence 'C') degraded i nto heat and ph otons 
via ion-neutral friction and cooling. iMullanl (jl97lh was 
first to discover that high magnetic fields can trans- 
fer kinetic energy to thermal energy in a c ontinu- 
ous manner and coined th e ter m C-shock. | Dr aine 
(|l980h . lDraine et al.l (|1983f h and lRoberge fc Drainel 
(1990) then described the multifluid nature of these 
shocks. C-shocks occur as long as the shock velocity 
remains below the propagation speed of the mag- 
netic precursor. This critical velocity above which a C- 
shock cannot exist is the magnetosonic velocity 



(8) 



where va c = B / y/Anp c is the Alfven speed in the charged 
fluid and c s is the speed of sound. The Alfven speed in 
the charges depends both on the magnetic field and on the 
inertia in the charged fluid. Thanks to the irradiation field, 
the gas is well ionised and provides ample free electrons to 
stick onto the grains. As a result most grains are charged 
negatively in our models and henc e they dominate the 
inertia in the charged fluid. Moreover. iGuillet et al.l l|2007f) 
showed that even the neutral grains spend an extensive 
fraction of their time attached to the magnetic field and 
they should also be included in the charged fluid inertia. 
We hence include all PAH and grains in the computation 
of the magnetosonic speed. As a result, the Alfven speed in 
the charged fluid vac is larger than va by a factor vac ~ 
yj p/pd va where p/pd ~ 180 is the gas to dust mass ratio. 
The actual number in our simulations turns out to be v m — 
21.2 km.s -1 , 21.7 km.s -1 and 22.9 km.s -1 for respectively 
nn = 10 2 cm -3 , 10 3 cm -3 and 10 4 cm -3 in the case 6=1 
and about ten times lower values for 6 = 0.1. 

Shocks with speed greater than v m will always be J- 
shocks. For instance, all shocks in our 6 = 0.1 grid are 
J-shocks. However, time-dependent shocks with velocities 
lower than v m could be C-shocks, J-shocks or even a com- 
binati on of the two (see IChieze et al.l [T998: Lesaffre " et al.l 
l2004al fq). Indeed, (i) the transverse magnetic fiel d could 
be lo wer for various orientations of the field (see IWardld 
Il998t for oblique C-shocks), (ii) or the C-shock could be 
at an earlier developme nt of its structure where it still 
has a J -shock component (jChieze et al.lll998l iLesaffre et al.l 
2004a), (iii) or the final fate of the shock could be a steady- 
state CJ-shock depending on the histor y of the shock 
(|Lesaffre et alJl2004aL ILesaffre et alJl2004b[) . As a result, we 
build our grids of models with C-shocks for shock velocities 
up to v m and we use J-shocks for higher shock velocities. 
In the following, we explore the properties of the resulting 
shocks. 



3.4. Weakly magnetized shocks (b=0.1) 

For high compression ratios (ie: Mach numbers much 
greater than 1), the maximum temperature in a J-shock 
(obtained just behind the leading adiabatic shock front) is 



2( 7 ~1) 



T - =11 tl-=-u 2 (9) 

^ ~ ( 7 + l) 2 k B [ > 

where p is the mean molecular weight of the gas, ks is the 
Boltzman constant and u is the shock velocity in the 
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shock frame. Note the shock velocity is very close to 
the velocity difference between upstream and down- 
stream gas for high Alfvenic Mach number. In par- 
ticular, this relation shows that the peak temperature in a 
shock is proportional to the square of the shock velocity. 

For such dilute molecular gas, only the lowest 
energy levels of H2 are populated and the appro- 
priate value for 7 is 7 = 5/3. With \i = 2.33 a.m.u. 
(which corresponds to the above = 1.4 in our near fully 
molecular conditions) relation (0) becomes: 

T max = 53K(u/lkm.s- 1 ) 2 . (10) 

Figure [IJa) displays the above theoretical value and the ac- 
tual maximum temperature reached in our models for each 
velocity tested in our grid at rtjj = 10 2 cm -3 . Both values 
remain close to one another. At high temperature (high 
velocity), the maximum temperature is closer to the value 
obtained in equation (O with 7 = 7/5. This stems from the 
excitation of rotational levels of the H2 molecules in the adi- 
abatic front. Note that we use a viscous length A = l/cmn 
with a viscous cross-section a = 3 x!0~ 15 cm 2 based on H2- 
H 2 collision data computed bv lMonchick fc Schaeferl (fl980) 
for a velocity dispersion of 1 km.s -1 . We checked that the 
viscous front between the peak temperature and the pre- 
shock is always adiabatic. Its width of 10 13 cm, visible 
on Figure [1(b), corresponds to the viscous length 
for 7Jh = 10 2 cm 3 . This figure surprisingly shows 
that higher velocity shocks return to pre-shock tem- 
perature on a smaller distance than low velocity 
shocks. However, the range of scales spans only one order 
of magnitude from 2 x 10 15 cm down to 2 x 10 14 cm with 
all shocks with velocities higher than 20 km.s -1 around the 
scale of 2 x 10 14 cm. Note that the range of flowing time 
through these shock structures also spans about one order 
of magnitude, from 400 yr to 2000 yr. For higher densi- 
ties, the picture is much the same and in particular time 
scales vary little within one category (J-type or C-type) 
of shocks. However, shocks which strongly dissociate H2 
(such as dense shocks at moderate velocity) have a wide 
H2 reformation plateau where H2 refo rmation provides th e 
necessary heat to keep the gas warm (jFlower et al.l l2~003). 

Figure shows the total column-density of H nuclei 
(Ah) across each shock in the grid. Also displayed is the 
total H2 column-density. The total A^h column-density is 
surprisingly constant around Ah=3-4 x 10 18 cm -2 over the 
range of velocities in our grid except at both ends. The in- 
crease of the compression factor at higher shock speed is 
compensated by the shorter cooling length. This result is 
in con trast to the results obtained by iHollenbach fc McKed 
(|1989f) in the velocity range 30 to 150 km.s -1 where the 
cooling column-density behind the shocks is seen to vary 
greatly. Shocks with velocities below about 7 km.s -1 have 
a lower than average total column-density, and shocks with 
velocities above 35 km.s -1 dissociate H2 which decreases 
the amount of cooling and increases the total column- 
density. 

Figures [UJb) and |2Jc) show the total column-densities 
across the shock for various neutrals and ions. Remember 
the end of each shock is defined when the temperature de- 
creases again under 20% above the pre-shock temperature. 
Since the overall column-density is a slowly varying func- 
tion of the shock velocity (figure [^a)) these plots also give 
a good estimate on how the average relative abundance of 




E_J , L_ ... 1. 

5 10 15 20 25 30 35 40 




10 12 10 13 i0 14 10 15 



distance (cm) 

(b) 

Figure 1. (a) Maximum temperature in the weakly magne- 
tized shocks, (b) Temperature profiles for some representa- 
tive weakly magnetized shocks, the fluid flows from left 
to right with the pre-shock on the left and the post- 
shock on the right. Only J-shocks are present for 
this low value of the magnetic field, nn = 10 2 cm -3 . 



each of these species varies with respect to the shock veloc- 
ity. One striking feature of these figures is the sharp 
rise of some molecular abundances at low veloc- 
ity. We trace these apparent velocity thresholds to 
temperature barriers or endothermicities of key re- 
actions which initiate the production of molecules. 
Equation (fTU|) allows to relate these critical temper- 
atures to key velocities in J-shocks. We list in table 
(|3|) the main bottleneck reactions which involve H2 
and a single atom or monoatomic ion along with 
their characteristic temperatures and J-shock ve- 
locities. At high velocities, the decrease in molecular con- 
tent (around 35 km.s -1 ) corresponds to the absence of H2 
molecules due to dissociation in the shock front. 

Figure^c) illustrates the variation of a few ions of inter- 
est. The abundance of C + is controlled mainly by photo- 
ionisation of C and charge exchange with PAHs and H2. 
As the velocity increases, the post-shock density increases 
with enhanced charge exchange, thus decreasing the C + 
abundance. For very large velocities (above 35 km.s -1 ) H2 
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Figure 2. For each shock velocity in our weakly magne- 
tized shocks (b = 0.1), we plot (a) the total H nuclei and 
H2 molecules column-densities (b) some neutrals of interest 
and (c) some ions of interest, njj = 10 2 cm' 3 . 



starts to be dissociated and C + increases again. In a molec- 
ular medium, the abundance of is controlled by the 
balance between cosmic ray ionisation of H2 (followed by a 
hydrogen exchange with H2) and recombination of H + . The 
charge balance links the ionisation degree to the abundance 



Reaction 


Temperature barrier 


Velocity 


0+H 2 - 


^OH+ H 


2980 K 


7.5 km.s -1 


C++ H 2 - 


->CH++ H 


4640 K 


9.4 km.s" 1 


S+ H 2 - 


->SH+ H 


9620 K 


13.5 km.s" 1 


S++ H 2 - 


->SH++ H 


9860 K 


13.6 km.s" 1 


C+ H 2 - 


->CH+ H 


14100 K 


16.3 km.s" 1 


Si++ H 2 - 


^SiH++ H 


14310 K 


16.4 km.s" 1 


N+ H 2 - 


->NH+ H 


14600 K 


16.6 km.s" 1 


H 2 dissociation energy 


52000 K 


31.3 km.s" 1 



Table 3. List of bottleneck reactions for molecules forma- 
tion with their temperature barrier or endothermicity and 
the corresponding velocity of the J-type shock that is able 
to provide such temperature computed using equation (|10j) . 
N + +H2 has little or no barrier: 168K. 



of C + , thus the abundance of H3 is inversely proportional 
to C + with an extra dependence on the square root of the 
temperature due to the temperature sensitivity of the re- 
combination rate. 

On figurc[3]wc display the variation of the total emission 
of various lines of interest depending on the shock velocities. 
This can be of interest for observers who want to charac- 
terize the velocities of observed shocks. Figure (3{a) shows 
some of the most observed atomic lines with the H2 0-0S(l) 
line for comparison. The C + ion emission decrease at mod- 
erate velocities is an abundance effect. Atomic O cooling is 
enhanced at moderate velocities because these shocks dis- 
sociate H2: this slows down the cooling and yields a tem- 
perature plateau where O is the dominant cooling agent. 
The stronger the shock, the more time it takes for H2 to 
reform, then the plateau widens and the O emission gets 
larger. 

Figure \3j(b) shows how the lowest energy lines of H2 
change with the velocity of the shock. The relative strengths 
of these lines vary greatly between u=3 km.s -1 and u=10 
km.s -1 . For example, the ratio S(l)/S(3) decreases from 
more than a hundred at 3 km.s -1 to about a fifth at 10 
km.s -1 . Indeed, the shock peak temperature in this range 
of velocities spans the energies of the upper levels of these 
transitions ( J = 7 corresponds to 3474 K, for instance, and 
see figure [lja) for how the maximum temperature in the 
shock varies with shock velocity): the excitation diagram 
in the low levels of H2 is greatly affected in this range of 
temperatures. By contrast, the aspect of the excitation di- 
agram does not change above 10 km.s" 1 : the emissivities 
of these lines are unable to discriminate the shock velocity 
u for u > 10 km.s" 1 . Indeed, above 10 km.s -1 the temper- 
ature experienced in these shocks is much greater than the 
J = 7 energy. Observers should use lines with upper levels 
with higher J values (such as rovibrational lines) to probe 
the velocities of these shocks or they should use diagnostics 
based on atomic lines. 

On figure U we examine what fraction of the kinetic 
energy flux input 

kinetic = \p u * ( n ) 

into the shock is radiated away by each coolant. Individual 
coolings are integrated through the whole shock structure 
as follows: 

hoo\= r Mz 
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Figure 3. (a) Atomic and (b) pure rotational H2 line emis- 
sion for weakly magnetized shocks, nn = 10 2 cm~ 3 . 



where A is the local rate of cooling and x\ is the end point 
of the shock (ie: where the temperature decreases back to 
20% above the pre-shock temperature). We also display 
the magnetic energy flux accross the shock 



h 



B flux 



Si, 

Air 



(13) 



where Ui is the ion velocity and B e is the magnetic 
held at the end of the shock. H2 cooling dominates 
almost everywhere except: (i) at very low velocity, where 
atomic O or C + cooling takes over depending on the density, 
and (ii) above the velocity for H2 dissociation (which is 
lower for higher densities) where atomic O and H (Lyman- 
a) cooling take over. At density njj = 10 4 cm~ 3 , and below 
20 km.s -1 , note that H2O cooling becomes important. 

3.5. Strongly magnetised shocks (b=l) 

In a C-shock, the kinetic energy is continuously transformed 
into thermal energy via ion-neutral friction. As a result, the 
heating is spread out on a much larger region than for J- 
shocks and the peak temperature is hence much smaller 
as seen on figure [SJa). However, reactions between neutral 
and ion species benefit from the ion- neutral drift, and the 



b=0.1,n M =10~ : main coolings 




Figure 4. Cooling integrated along the shock length, nor- 
malised by kinetic power \pv? in weakly magnetized shocks, 
for nn = 10 2 cm~ 3 . The solid line with triangles shows the 
fraction of the power which is transferred into a flux of 
magnetic energy. On this figure, H2 cooling incorporates 
the cooling from 200 lines between the 149 H2 levels in- 
cluded in the simulation, whereas figure [3] considers only 6 
individual lines amongst the lowest levels of H2 . 



effective temperature which is used to compute the reac- 
tion rate is higher (see iPineau des Forets et alj |l986). For 
instance, the effective temperature for reactions involving 
C + and H2 is 



T eff (C+, H 2 ) 



2T„ + 12T, 2 x 12 m r 



2 + 12 



2 + 12 3k B 



{u n - Ui f (14) 



where (u„ — u{) is the local drift velocity between the ions 
and the neutrals, m p is the proton mass and T n and Tj are 
respectively the temperatures of the neutral and ionised 
gas. Figure 0(a) shows the maximum effective temperature 
of the reaction C + + H2, which incidentally is very close 
to the maximum temperature which would be obtained in 
the corresponding J-shock. 

Compared to the temperature profile of a J-shock, the 
temperature profile of a C-shock is broad, with a lower 
temperature: see figure E^b). However, within each type 
of shock, the total size of the shock structure varies over 
less than an order of magnitude for the range of velocities 
we tested, as is apparent on figure [SJb). This is even more 
striking for the flowing time scales across the shock struc- 
tures which, for tih = 10 2 cm" 3 , varies only from 8000 to 
6000 years for the C-shocks and from 1500 to 1000 years 
for the J-shocks. As mentioned above, this holds for higher 
densities as well. 

Slightly above 20 km.s" 1 , all shocks in our grids of mod- 
els are J-shocks, but the strength of the magnetic field varies 
from b = 0.1 in the previous subsection to b = 1 in this 
subsection. The higher magnetic field limits the compres- 
sion in the shock and the collisional processes take longer 
to occur. In particular, the cooling length for the J-shocks 
above 20 km.s -1 is between 4 x 10 15 cm and 7 x 10 15 cm: 
much wider than for the corresponding weakly magnetised 
J-shocks. The lower density but larger size of the shock im- 
pacts on the chemical composition and structure in subtle 
ways. A careful comparison of the right hand sides of figures 
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Figure 5. (a) Maximum temperature in the highly magne- 
tized shocks (b) Temperature profiles for some representa- 
tive highly magnetized shocks, the fluid flows from left 
to right with the pre-shock on the left and the post- 
shock on the right, njr = 10 2 cm -3 . 



[2] and |5] shows the column-density of H nuclei is higher in 
magnetised shocks and moderate variations in the chemical 
composition are noticeable. 

Figure [5] shows the molecule production over the range 
of models in our grid at njj = 10 2 cm -3 . The transition 
from C-shocks (velocity lower than 21 km.s -1 ) to J-shocks 
(higher shock speed) is clear. The column-densities of most 
species have a drop at 21 km.s -1 which expresses a boosted 
molecular production in C-shocks. Three effects are at play 
which favour molecules production in C-shocks. 

First, the ion-neutral drift helps to overcome reaction 
barriers of ion-neutral reactions (which impact C-bearing 
species, for example, whose bottleneck reaction is C + +H2). 
Second, the resulting frictional heating keeps the tempera- 
ture warm throughout the shock. By comparison a J-shock 
is very hot at the peak temperature right behind the vis- 
cous front but quickly cools down in the trailing relaxation 
layer. Third, the steady frictional heating slows down cool- 
ing and compression in the relaxation layer and the net 
result is an increased total Nn column-density in C-shocks 
compared to J-shocks (see the drop at 21 km.s -1 in figure 
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Figure 6. Same as figure [5] for highly magnetized shocks 
(b = 1). n H = 10 2 cm' 3 . 



[6^,). This is the main factor which impacts those molecular 
species whose production relies mainly on neutral-neutral 
reactions (such as O-bearing species). In particular, these 
species show a more gradual rise at low velocities, because 
of the slower rise of the maximum neutral temperature in 
C-shocks. 
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Figure 7. (a) Atomic and (b) rotational H2 lines emission 
for highly magnetized shocks, nn = 10 2 cm -3 . 



We display the total integrated emission in various lines 
of interest on figures [TJa) and(7Jb): as for figure [31 these 
give clues on which observations constrain what shock ve- 
locity. In particular, the emissivities of H2 lines with low J 
upper levels now vary from 3 to 20 km.s , but are indepen- 
dent of velocity above 20 km.s -1 . The neutral temperature 
in C-shocks rises more slowly than J-shocks as shown by 
the comparison between figures [T] and [SJ but the maximum 
temperature in C-shocks is more representative of the tem- 
perature in the whole mass of the shock and there is more 
column-density in C-shocks. As a result, the left hand sides 
(between 3 and 20 km.s -1 ) of figuresJTJa) andJTJb) look like 
blow ups of figures Eta) and EJb) shifted to greater emis- 
sivities. However, C + stands out even more for C-shocks 
because its excitation benefits from the ion- neutral drift. 
The right hand sides of figures [TJa) and [TJb) are similar 
to the corresponding weakly magnetised J-shocks seen in 
figures G3a) aiid[3]Jb): the emissivities are poorly affected 
by the density change due to the magnetic field except for 
C+. 

In figure [5] we examine, averaged over the whole shock 
structure as in figure [4J what fraction of the kinetic energy 
flux input into the shock is radiated away by each coolant 
in the low density case (rin = 10 2 cm -3 ). At very low ve- 
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Figure 8. Cooling normalised by kinetic power in highly 
magnetized shocks for pre-shock density njj = 10 2 cm -3 . 
The solid line with triangles shows the fraction of the power 
which is transferred into a flux of magnetic energy. 

locity, C + is the most efficient coolant (rather than O for 
weakly magnetised shocks), but H2 takes over for shock 
velocities shortly above 5 km.s -1 . Molecules such as CO, 
H2O or OH carry less than a percent of this energy flux. 
The kinetic energy flux is mainly transferred into a mag- 
netic energy flux (via field compression) for velocities below 
10 km.s -1 and is mainly radiated away above that veloc- 
ity. At higher densities, the low velocity coolant becomes 
O. It is also interesting to note that at nn = 10 4 cm -3 
H2O becomes an important cooling agent for low velocities. 
iFlower fc Pineau des Foretsl (|2010h also quoted the fraction 
of energy radiated away in various cooling processes in their 
shock models. However, in the present work the irradiation 
field photo-dissociates CO and ionises C so that C + cooling 
is more prominent and CO cooling is less important than 
in their work, but the results are otherwise similar. In par- 
ticular, they also find that the magnetic energy flux stores 
50% of the injected power at u — 10 km.s -1 . 

4. PDF of shocks 

4.1. Method of fit 

Numerical simulations of driven supersonic turbulence show 
that the med i um ex periences a whole range of shock speeds: 
ISmith et~ al. (2000) show that the probability distribution 
function (PDF) of the velocity jumps follows a decreasing 
power law with an exponential cut off at a Mach number of 
a few. iPetv fc Falgaronel (|2000h carefully separated vortical 
from compressional contributions in the velocity field of a 
compressible simulation of decaying supersonic turbulence: 
they uncover an exponential distribution for the conver- 
gence time scales — l/div(w) where v is the fluid velocity. 
Subsonic wind tunnel experiments display an expone ntial 
distribution of velocity increments (jMouri et al.l l2008). 

These results hint at a variety of possible PDF shapes 
for the statistics of s hocks, all favoura b ly wei ghted towards 
low velocity shocks. iGuillard et al.l ((2009) have mod- 
eled Spitzer H2 observations of SQ with only two 
discrete values of the shock velocities. Here, we aim 
at fitting a statistical distribution of shocks to a collection 
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of observed quantities. Guided by the above remarks on 
turbulence, we fit various shapes of shocks PDF: decaying 
exponentials, power-laws, Gaussians and a piece- wise expo- 
nential fit. 

We consider among j — 1..M observable quantities 
an observable Xj{u) associated to a steady shock of ve- 
locity u (expressed in km.s -1 ). We note f(u, {pi}i=i..N) 
the probability density parametrised by the set of parame- 
ters {pi}i = i..N ■ Of course, the fitting makes sense only as 
long as the number of observed quantities M is sufficiently 
greater than the number of fitted parameters N. For exam- 
ple, an exponential PDF reads 

fe(u,Pi) = C(px) exp(-piu) (15) 
with the normalisation constant 

C(pi) = l/Y, exp(-pm). (16) 

We then compute observable quantities averaged over 
the PDF of shocks as 

In order to ease the fitting process, we normalise the 
above values to the sum of the observable quantities: 

V, = ^- M f X 3 (18) 

where {X°}j = i..M are the actual observed values. Note that 
this normalisation renders irrelevant the normalisation con- 
stant of each PDF, but it decreases by one unit the number 
of degrees of freedom Aff rco = M — N — 1. We finally find 
the optimal set of parameters {p^},-i..Ar which minimises 
the quantity 

X 2 ({kW..aO = j^— E (logio^-logxo^ ) 2 /^ 2 
Ju free j=l „ M 

where Wj is the uncertainty on the observed value log 10 X®. 
We also retrieve estimates {£'i}i = i..jv of the errors on the 
parameters by computing the diagonal elements of the co- 

variance matrix of -§^iX 011 t ne optimal set. 

Table 2] summarises the various PDFs we tried and their 
parametrization. The exponential and power-law PDFs are 
motivated by the above mentioned work on the statisti- 
cal distribution of shocks in numerical and laboratory ex- 
periments. The piece-wise exponential PDF is built on in- 
tervals designed such that the observable quantities vary 
significantly (cf. figures [3] and [7]) . Otherwise, the parame- 
ters would be degenerate. We use as parameters the PDF 
value at the transition between these intervals. Thus, these 
parameters' optimal values and corresponding error bars 
directly probe the PDF and its uncertainty at these chosen 
locations. Finally, we use a single or a double Gaussian PDF 
to mimic the fit of one shock or two independent shocks: th is 
allows to connect with the work of iGuillard et al.l (|2009f ). 

In addition, we note that the choice of the boundaries 
(upper and lower velocities) for the range of our grid of 
models does not affect significantly our results. This stems 



PDF formula N 

Power-Law u~ pi 1 

Exponential exp(— p\u) 1 

„. . at u = 3, 10, 20, 40 

Piece-wise exponential cl i 3 

/ W = 1, Pi, f>2, P3 

1-Gaussian e -(«-Pi) 2 i 



2- Gaussian e~ ( "~ pi)2 + p 3 e~ ( "~ ra)2 3 

Table 4. Template PDFs adjusted to the observed data 
(normalisation constants are discarded as irrelevant, see 
text). 



from the fact that the shock emission is very small at low 
velocity and high velocity shocks are very rare in the so- 
lutions we obtain. Finally, grids of models with velocity 
spacings of 1 km.s^ 1 and 0.5 km.s -1 gave similar results, 
which validates the resolution we chose. 

4.2. Applications 

4.2.1. The Stephan's Quintet galaxy collision 

In this section we use our model results to fit H2 line emis- 
sion from the galaxy-wide shock in SQ. This shock struc- 
ture, first identified in radio continuum observations and 
X-ray images was discovered to be a luminous H 2 sou rce 
with Spitzer (|Appleton et al.l 120061 : ICluver et al.ll2010D . It 
is associated with the entry of a galaxy into a group of in- 
teracting galaxies with a relative velocity ~ 1000 km.s -1 . 
The luminosity of the molecular gas in H2 rotational lines 
is observed to be larg er than that of the plasma in X-rays. 
IGuillard et alJ (|2009T ) present a first interpretation where 
the H2 emission is powered by dissipation of turbulence 
driven by the large scale collision. CO observations have 
since shown that the kinetic energy of the molecular gas is 
larger than the thermal energy of the hot (X-ray emitting) 
plasma. It is the mai n energy reservoir av a ilable to power 
the H 2 line emission (|Guillard et alJ I2012T ) . IGuillard et al.l 
(2009]) show that the H2 excitation is well fit by a combina- 
tion of two MHD shocks with velocities of 5 and 20km.s~ 1 
in de nse gas (ny\ ~ 10 3 to 10 4 cm~ 3 ), us ing the shock mod- 
els of lFlower fc Pineau des Foretsl (|2003f ) , which do not take 
into account UV radiation. Our grid of models allows us to 
test an alternative interpretation where the H2 emission is 
accounted for with MHD shocks in lower density UV irra- 
diated gas. In doing so we quantify a solution where the 
physical state of the H2 gas in SQ is akin to that of the 
cold neutral medium in the Galaxy. Evidence that the gas is 
magnetized comes from radio continuum observations. The 
synchrotron brightness of the shock yields a magnetic field 
value B ~ 10 p,G, assuming equipartition of magnetic and 
cosmic-ray energy (Xu et al. 2003) . This value is compara - 
ble to the field strength reported bv lCrutcher et al. (2010) 
for the Galactic diffuse ISM. The mean UV radiation field 
in th e SQ shock estimated to be Go = 1.4 (jGuillard et alJ 
2010), is also close to the reference value for the ISM in the 
Solar Neighborhood. 

The wide extension (30 kpc) of the SQ shock suggests 
that in the relaxation layer which follows the main high 
velocity (~ 600 km.s" 1 ) shock, a large number of sub- 
structures are formed which then collide and yield a range of 
sub-s hocks with much smaller velocities (see IGuillard et all 
2009). Simulations of large scale shocks subject to ther- 
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mal instability indeed show that turbulence is sustaine d in 
the t railing relaxation layer (Kovama fc Inuts uka 200 j and 
later I Audit fc Hennebelldl2005t Irlennebelle fc Audit! l2007t 
lAudit fe Hennebelld l2010f) . This connects with the theo- 
retical work on turbulence mentioned above, and justifies 
our investigation with continuous PDFs of shock velocities. 
The observable targets which we aim at reproducing are 
the emiss ion values for the si x H2 lines 0-0 5(0) to 5(5) in 
table 1 of |a uver et all (|2010D which are integrated over the 
main shock structure. 

Table [5] shows the best \ 2 obtained in all our attempts 
to fit the data. Two solutions stand out with \ 2 values rea- 
sonably close to one, namely: the piece-wise exponential 
and the two Gaussians at b = 1 and nn = 10 2 cm -3 (their 
X 2 values are emphasised in bold faces in table [5]). Figure 
[9] illustrates the quality of the fit on the observed H2-lines. 
Our models reproduce quite accurately every H2-line ex- 
cept for 0-0 5(4)Q. For comparison, when we use th e same 
technique and same grid of models as lGuillard et al.l (|2009t ) 
(ie: models without irradiation at nn = 10 4 cm -3 ), we find 
the best fit is obtained for two shocks of velocity 5 and 22 
km.s -1 in proportion 1:0.008 with a reduced \ 2 = 78 (com- 
pared to x 2 = 15-8 for velocities 3.7 and 21 km.s -1 in pro- 
portion 1:0.0002 in our 2-Gaussians high-density case). The 
coincidence of the optimal parameters with our 2-Gaussian 
solution is surprising, but the improvement on the \ 2 m 
the present work is also striking. Irradiation does really im- 
prove the comparison with the observations. Our work also 
shows that low-density solutions are even more viable than 
the previously found high-density solutions. 

All optimal PDF shapes we find (including those with 
a large % 2 ) are statistically biased towards low velocity 
shocks. Figure [TD] displays the two best fit PDF solutions. 
The error bars (gray regions in these figures) are deter- 
mined as follows: we vary the parameters of the PDFs at 
the extreme of their 3-cr range of uncertainty and we take 
the minimum and maximum values predicted when using 
these extreme PDFs. The piece-wise exponential adjust- 
ment shows a dip at 10 km.s -1 which indicates a bimodal 
distribution of shocks consistent with the two-Gaussian fit 
(which is bimodal in essence...). Perturbation tests of this 
fit show that the PDF level at 10 km.s -1 depends mainly 
on the ratio S(0)/S(1), with a deeper dip for higher val- 
ues of this ratio. On the other hand, the currently observed 
data do not discriminate whether the distribution of mod- 
erate velocity shocks is wide spread or well centred on a 
distinct velocity. Indeed, as seen on figure [7Jb), above 20 
km.s -1 the emission properties of the observed H2 lines do 
not change with velocity, hence these lines cannot probe the 
shape of the PDF in this range of velocities. Ro- vibrational 
lines should be used to probe the shocks at these veloci- 
ties. We note all piecewise-exponential solutions (including 
those with a larger \ 2 & t higher densities) are similarly 
consistent with bimodal distributions. This is not expected 
from the statistics of shock velocities in numerical simu- 
lations of turbulence. It may be a result of the intrinsic 
multiphase nature of the interstellar medium, where the 
moderate velocity shocks would be interpreted as collisions 
between molecular clouds and low velocity shocks would be 
the signature of turbulence dissipation within these clouds. 



2 This line actually lies at the position of the 7.7pm PAH 
emission feature, which makes it difficult to measure accurately. 



For each of the optimal PDF found, we can predict the 
value X of other quantities of interest. We provide in ta- 
ble [5] the total column densities for some molecules and the 
expected emission of some atomic lines in our best two so- 
lutions. These total column-densities are integrated 
from the pre-shock to the point where the temper- 
ature decreases back to 20% above the pre-shock 
temperature. In particular, we provide a measure of the 
total H2 column-density: the piece-wise exponential model 
yields 

^shocked(H 2 ) = 2.6 x 10 20 cm -2 
and the 2-Gaussian model yields 

^shocked^) = 1-0 >< 10 2 °cm- 2 . 

Recent mea s ureme nts from CO spectroscopy 
(|Guillard et al.l |2012( ) allow to estimate the total H2 
column-density in the SQ large scale shock structure as 

iV total (H 2 ) = 8.5 x 10 20 ± 1.5 x 10 20 cuT 2 

(assuming a Galactic value for the CO emission to H2 col- 
umn conversion factor). This puts the fraction of shocked 
gas in this line of sight between 32% and 12%, depend- 
ing on which model w e adopt, which confirms the results of 
Guill ard et alJ (|2010f ) that quite substantial amounts of gas 
are shock-heated. In fact, quite a few of our optimal solu- 
tions predict A/ g j 10c T ie( j(H2) > 10 21 cm -2 and are therefore 
probably ruled out (only the viable solutions are underlined 
in table [5]). 

Since the shocks account for a significant fraction of the 
gas mass, any chemical enhancements in the shocks may 
have a significant contribution to molecular abundances. 
Table [6] also gives estimates of the chemical yields in these 
shocks. For example, the H2O abundance is on the order 
of 10 -6 , significantly greate r to that observed in the dif- 
fuse ISM in the Milky- Way (|Wvrowski et al.ll2010h . These 
abundances could be even higher if the medium is clumped 
(and therefore shielded) and Herschel may provide a good 
opportunity to test this. For standard irradiation C + is the 
dominant carbon gas species and this is also a result which 
will depend on shielding. 

The importance of these species (in addition to H2), 
as cooling agents of the interstellar turbulence depends di- 
rectly on their abundance. Herschel observations could thus 
provide additional constraints to achieve a precise mod- 
elling of the physical and chemical state of the gas. We 
provide estimates of the total integrated flux in the molec- 
ular cooling agents in table [5] H2 itself is responsible for 
more than 94% of the cooling in both our best-fit mod- 
els. However its emission is distributed over many different 
lines, the strongest of which H2-S(l) contributes by a few 
percent of the total. Note that for low-velocity C-shocks 
the pure rotational lines sum up to almost all the H2 cool- 
ing. On the contrary, moderate velocity J-shocks experience 
much higher temperatures and the energy in these shocks is 
radiated away through higher H2-levels. If our models with 
a contribution from moderate velocity J-shocks hold, a lot 
more energy should be emitted from higher excitation levels 
of the H2 molecule than is observed from pure rotational 
lines. Finally, note that other molecular coolants contribute 
at most a few percent of the H2-S(l) luminosity. 

The emission of the C + line at 158 /jm is quite strong 
in our best two models (about half as strong as the H2-S(l) 
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b 


n H 


1-Gauss 


pow-law 


exp. 


pw-exp. 


2-Gauss 


0.1 


W 2 


371.8 


2307.0 


54.3 


60.8 


11.2 


0.1 


10 3 


504.0 


1650.4 


152.4 


61.1 


105.6 


0.1 


10 4 


416.1 


2139.9 


174.3 


580.8 


155.3 


1 


10 2 


1628.5 


184.2 


598.5 


2.6 


2.0 


1 


10 3 


139.3 


175.1 


35.9 


5.0 


13.8 


1 


10 4 


130.3 


1648.0 


12.6 


6.3 


15.8 



Table 5. Summary of the optimal x values obtained 
in all our attempts to fit the SQ data. Both our opti- 
mal values are emphasised in bold faces, hh is given in 
cm~ 3 . Underlined values correspond to solutions which 
have A(H 2 ) < 



H 2 rotational lines: b=l,n u =W" 




line), and so is the OI emission to a lesser extent. Assuming 
all carbon is in the form of C + in the un-shocked fraction of 
the total line of sight column-density N^^i^) , we pre- 
dict that the C + line should shine about as much as the H2- 
S(l) line, with a significant contribution from shocks (66% 
and 52% respectively in both models). Interestingly, this 
was also sugges ted in an inte r preta tion of AKARI 
observations by I Suzuki et ahl (12011 b In contrast, the 
model proposed by iGuillard et all (120091) with no UV ir- 
radiation has no CII emission. This means that C + can 
indeed be a good signature for the dissipation of kinetic 
and magn etic energy in weakly shielded gas (as previously 
found by iFalearone et al.ll2007h . We also note that our vi- 
able models at higher densities predict even stronger emis- 
sion in C + by up to a factor of three. As such, C + emis- 
sion measurements will help to probe both the density and 
shielding of the gas. By contrast to C + , the measured 
emission of the Si + line at 34.8 /im is much greater 
than the predicted emission from our shock mod- 
els. It is probably dominated by the contribution 
fro m the hot ion i sed m edium, as already suggested 
bv lCluver et alJ (|2010h . 

Figure [11] displays the distribution of energy radiated 
away for several cooling agents in the best-fit 2-Gaussian so- 
lution. A significant fraction of energy is dissipated in low- 
velocity shocks but more energy is radiated within moder- 
ate velocity shocks. Indeed the dissipated power goes as the 
cube of the velocity, which slightly more than compensates 
the higher number of low- velocity shocks. The figure clearly 
demonstrates that H2 is by far the main cooling agent. It 
also shows that the emission from atomic cooling agents 
comes mainly from low- velocity shocks whereas the molec- 
ular emission is dominated by moderate velocity shocks. 
Similarly, we have checked that the emission from low J 
H2-lines probe low-velocity shocks whereas higher J lines 
probe larger velocity shocks, which seems rather natural. 

4.2.2. Chamaeleon 
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iNehme et al.1 (|2008l ) successfully interpreted observational 
results of a line of sight sampling diffuse molecular gas 
in front of the star HD 102065 in Chamaeleon. Their 
PDR model with G = 0.4, n H = 80 cm~ 3 and A v = 0.67 
reproduced accurately as many as seven independent obser- 
vational quantities. However, they were unable to account 
for the column-densities of rotationaly excited states of H2 
and failed to reproduce the observed abundance of CH + . 
Here, we assume that both a PDR and a statistical distri- 
bution of shocks contribute to this line of sight. Hence we 



Figure 9. Observed f luxes of the H2 lin es in the SQ (red 
dots with error bars: [a uver et alJ I2010L see also table EJ} 
with the results of our best two models. Note most error 
bars are so small they remain barely visible. 



PDFs for 6=1, n„ =10 J 



2-Ganss1aTi 




15 20 25 

u (km/a) 

Figure 10. Shape of our best two optimal PDFs fit to the 
SQ data. Grey areas represent uncertainties due to the 
propagation of observational errors (see details in text). 



attempt to fit as many as 11 observational measurements 
available for this line of sight with a linear combination of 
our shock models and a PDR model. 

The PDR model is computed with our shock code in 
a PDR mode, see section 12.11 In principle, we should fit 
all irradiation parameters of the PDR alongside the grid of 
shocks. But this would assume we recompute the whole grid 
of shocks each time we probe new irradiation parameters 
and would lead to prohibitive computational time. Hence 
we fixed the irradiation conditions and the PDR model 
therefore introduces only one extra parameter in the fit: 
its weight in the line of sight compared to the shock mod- 
els. We use Go — 4 and A v — 0.335 as entrance conditions 
for both the PDR and the shocks and we compute a new 
grid of shock models from 3 to 40 km.s -1 with nn = 80 
cm~ 3 and an ortho/para ratio of 0.7 as is observed. 

We thus obtain reasonable x 2 values with a best value of 
11 for both the 2-Gaussian and the piecewise-exponential. 
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name pw-exp 2-Gauss obs. err. 



H2-lines observations (erg. cm 


*^.s ^ .sr 


) 


0-0S(0) 


5.7(-7) 


5.6(-7) 


5.8(-7) 


3.5(-8) 


O-OS(l) 


4.2(-6) 


4.2(-6) 


4.3(-6) 


4.8(-8) 


0-0S(2) 


1.7(-6) 


1.7(-6) 


1.7(-6) 


7.0(-8) 


0-0S(3) 


3.8(-6) 


4.0(-6) 


4.2(-6) 


1.5(-7) 


0-0S(4) 


9.0(-7) 


9.4(-7) 


4.7(-7) 


1.8(-7) 


0-0S(5) 


2.5(-6) 


2.4(-6) 


2.6(-6) 


1.3(-7) 


Predictions 


Column-densities (cm -2 ) 


H 2 (0,0) 


5.8(19) 


2.2(19) 


- 


- 


Ha(0,l) 


1.8(20) 


7.6(19) 


- 


- 


H 2 (0,2) 


3.4(18) 


3.4(18) 


- 


- 


N(H a ) 


2.5(20) 


1.0(20) 


<1.0(21) 


- 


N(H) 


4.3(19) 


1.8(19) 


- 


- 


N(C+) 


7.6(16) 


3.2(16) 


- 


- 


N(CO) 


2.5(14) 


1.8(14) 


- 


- 


N(H 2 0) 


1.6(14) 


1.9(14) 


- 


- 


N(OH) 


1.9(14) 


1.7(14) 


- 


- 


N(CH) 


2.1(13) 


1.0(13) 






N(CH+) 


1.7(12) 


1.1(12) 






N(HCO+) 


7.4(11) 


6.2(11) 






Atomic lines (erg. cm z .s 


-^sr" 1 ) 




S(25.2m) 


1.4(-8) 


7.0(-9) 






Si+(34.8m) 


9.8(-8) 


7.9(-8) 


1.6(-6) 


5(-8) 


0(63.2m) 


1.4(-6) 


9.5(-7) 






0(145.3m) 


1.2(-7) 


9.0(-8) 






C+(158m) 


3.4(-6) 


2.2(-6) 






C(370.4m) 


4.6(-8) 


3.1(-8) 






C(609.8m) 


2.2(-8) 


1.3(-8) 






Integrated molecular cooling (erg 


;.cm _;! .s" 




E(H 2 ) 


1.6(-4) 


5.4(-5) 






E(CO) 


4.4(-8) 


2.9(-8) 






E(H 2 0) 


l.l(-7) 


1.0(-7) 






E(OH) 


4.1(-7) 


1.0(-7) 







Table 6. Results of our best two fits to the observed 
spectral line energy distribution in the SQ (flux are given 
in erg.cm~ 2 .s .sr -1 ) and their respective predictions: 
column-densities are in cm -2 and cooling is indicated in 
erg.cm .s .sr . Parentheses denote powers of ten. 



The error mainly comes from overestimating the column- 
density in J = 5, but the molecular chemistry of the line of 
sight is reproduced rather accurately (see table [7J . 

The best two fitting PDFs compare qualitatively well 
to the solutions we found for SQ, except the low veloc- 
ity part of the 2-Gaussian is more tightly peaked on low 
velocities (the moderate velocity component accounts for 
3% of the total probability), and the piecewise-exponential 
dip at 10 km.s" 1 is less pronounced. The two PDF shapes 
are hence slightly less consistent than for the SQ, and con- 
tinuous solutions are not ruled out (the exponential and 
power-law solutions yield x 2 values on the order of 20 and 
15 respectively). However, in the case of the Chamaeleon, 
the physical extent of the measured beam is much smaller 
than for the SQ. In particular, the total number of shocks 
involved is presumably much smaller and discrete numbers 
effects might show up in the distribution. 

Our best two models predict similar values for most of 
the observables listed in table [7J We now have two contri- 
butions from the PDR and the shocks and for each quantity 
listed in table [JJ we quote in brackets which fraction comes 
from the shocks. The weight of the PDR model in both 
best solutions is slightly different, and this yields a mass- 




Figure 11. Distribution of the energy dissipated via several 
cooling agents for the 2-Gaussian best-fit model. 



fraction of shocked matter on the line of sight of 4% in the 
piecewise-exponential model and 19% in the 2-Gaussians 
model. The contribution from shocks is hence small but 
not insignificant. However, molecules other than H2 are all 
formed in the shocks and not in the PDR, with for example 
98% of the CO coming from the shocks. Note that in both 
our two best models, shocks with velocity above 15 km.s -1 
account for the formation of almost all molecules. This is 
due to the overcome of chemical barriers thanks to the ther- 
mal energy released in shocks and to the ion-neutrals drift, 
as mentionned in sections 13.41 and 13.51 . 

CH + needs mechanical energy injection to raise the ef- 
fective temperature for its formation which is why pure 
PDR models failed to account for its abundance. In fact, 
CH + forces a too hot component into the model, which 
then over-produces the J = 5 H2-level. A possible solution 
would be that the irradiation in the far-UV is enhanced in 
the Chamaeleon compared to the Draine ISRF: this would 
make more CH + from CH^ and CH^" photo-dissociation 
(jFalgarone et al.ll2010h while requiring a lower temperature 
more compatible with the obs erved J = 5 H2-level col umn- 
density. In fact, UV spectra in lBoulaneer et all (|l994l) sug- 
gest that the enhanced far-UV field could be accounted by 
a flatter extinction curve in the far-UV for this very line of 
sight. 

As for the SQ, H2 is by far the main cooling agent (it 
radiates more than 92% of the total cooling in both best-fit 
models) but CII (and OI to a lesser extent) has emission 
comparable to the H2-S(l) line. Even though it is now dom- 
inated by the emission from the PDR, the shock contribu- 
tion to CII emission is still significant (from 18% to 31% 
according to both best models). Interestingly, we predict 
OH should shine about as much as CII and as for the SQ 
solution, the OH emission is almost completely dominated 
by the contribution from shocks at moderate velocity. OH 
observations would hence provide a good test of the exis- 
tence or not of a moderate velocity shock component in the 
Chamaeleon. 
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piecewise-exp 2- Gaussian 



obs 



Column-densities (cm ) 



Observations 



"SQ shock). We fit the observables with a continuous com- 
bination of shock models as a phenomenological descrip- 



N(H) 
H a (0,0) 
H 2 (0,1) 
H 2 (0,2) 
H a (0,3) 
H 2 (0,4) 
H 2 (0,5) 

N(C) 
N(CO) 
N(CH) 
N(CH+) 



2.6(20 
1.8(20 
2.1(20 
2.6(18 
4.2(17 
4.4(16 
2.2(16 
2.4(15 
5.2(13 
6.2(12 
8.8(12 



N(H 2 ) 
N(C+) 
N(H 2 0) 
N(OH) 
N(HCO+) 
N(CN) 
N(C a ) 



3.9(20 
1.4(17 
2.8(13 
6.9(13 
2.7(11 
3.9(11 
3.6(10 



Atomic lines (erg 



S(25.2m) 
Si+(34.8m) 

0(63.2m) 
0(145.3m) 
C+(158m) 
C(370.4m) 
C(609.8m) 



4.5(-9 
6.4(-8 
1.2(-6 
8.7(-8 
3.8(-6 
2.2(-8 
1.3(-8 



Molecular coolint 



0-0S(l) 
E(H 2 ) 
E(CO) 
E(H 2 0) 
E(OH) 



1.9(-6 
9.5(-5 
9.3(-8 
1.8(-7 
2.9(-6 



O06] 
0.04] 
0.03] 
0.53] 
0.99] 
1.00] 
1.00] 
0.20] 
0.98] 
0.87] 
1.00] 



2.7 
1.8 
1.9 
2.6 
4.9 
5.9 
3.8 
2.3 
4.9 
4.8 
1.0 



Predictions 



0M[ 

0.05] 

0.99] 

0.99] 

0.99] 

0.99] 

0.80] 



3.8 
1.4 
5.3 
1.1 
3.4 
5.3 
2.8 



.sr 



0.82] 
0.66] 
0.76] 
0.18] 
0.60] 
0.33] 



5.1 
6.0 
1.1 
7.7 
3.9 
2.3 
1.3 



(erg. cm .s" 



099] 
1.00] 
1.00] 
1.00] 
1.00] 



2.2 
5.7 
8.5 
2.3 
2.0 



023] 
0.20] 
0.14] 
0.61] 
0.99] 
00] 
00] 
0.30] 
0.98] 
0.86] 
1.00] 



3.1(20) 
2.0(20) 
1.4(20) 
2.6(18) 
2.0(17) 
2.0(16) 
7.1(14) 
6.0(14) 
5.6(13) 
6.4(12) 
1.2(13) 



OT7T 

0.19] 

1.00] 

0.99] 

0.99] 

0.99] 

0.77] 



<5.9(11) 
<1.4(13) 



0.84] 
0.67] 
0.77] 
0.31] 
0.68] 
0.44] 



2.8(-6) 



099] 
0.99] 
1.00] 
1.00] 
1.00] 



Table 7. Results and predictions of the fit of piecewise- 
exponential to th e Chamaeleon data (|Grv et al.l [2002; 
iNehme et aD l2008h . Powers of ten are indicated in 
parentheses. For each value, we provide the fraction com- 
ing from the shocks in brackets. 



5. Summary, conclusions, discussion, prospects 

Weprovide to the scientific community a grid of shock mod- 
els in UV heated gas which may be used to interpret obser- 
vations of the Milky- Way diffuse ISM and integrated prop- 
erties of galaxies in general. We examine magnetised shocks 
in media with densities from njj = 10 2 to 10 4 cm~ 3 , with 
standard ISM irradiation conditions and at low to moder- 
ate velocities (from 3 to 40 km.s -1 ). When the velocity of 
these shocks is below the critical velocity for the existence 
of pure C-type shocks, we compute C-type shocks, other- 
wise, we compute J-type shocks. The models include the 
effects of the ambient UV irradiation field on the pre-shock 
chemistry and thereby on the relative importance of cooling 
lines in the shock. For instance, C + can provide significant 
line cooling in shocks propagating in the UV irradiated gas 
where it is the dominant carbon species. 

We illustrate how the model results may be used to 
interpret data on one galactic line of sight through the dif- 
fuse ISM (Chamaeleon) and one extra-galactic object (the 



The output data of our models is archived on 



http://cemag.ens.fr 



tio n of the complex statistical properties of the dissipation 
3-0(19jbf turbulence. The model provides a rather good match 
2-0(19j;o the data for bimodal distributions of shock velocities. 
1.0(19)r;iri s d oes not match pre dictions from numerical simula- 
, „>,„tions of lSmith et ail (|2000h who found a power-law distribu- 
1 0(16/* 011 w ^ n an exponential cut-off. Further work is needed to 
5 o(14 un< ^ ers ^ an< ^ t ne reasons behind this apparent discrepancy. 
1.5(14j- n terestingly, the low and moderate velocity components of 
6.8(12*pur best-fit PDF operate in very different regimes of energy 
1.0(12)iissipation. Indeed, in our best-fit models to the observa- 
1.9(12);io ns, the low velocity shocks are of C-type and they dissi- 
p_ate energy via ion-neutral drift whereas the larger velocity 

- shocks are of J-type and they undergo viscous dissipation. 

In both our interpretations of the SQ and the 
Chamaeleon, a significant fraction of the molecular gas 
is shock heated and the chemistry in this shock heated 
gas has a dominant contribution to key molecules such as 
CO, H2O, OH or CH + which are commonly used as diag- 
nostics. This shock heated contribution might be greater 
I th an the contribution of C O computed in simulations 
_ by iGlover fc Mac Lowl (|2011| ) where the resolution in the 

- shocks is too scarce to trigger the shock chemistry: too 

- low resolution smears out the temperature to such 

- an extent that shocks become nearly isothermal 
" with a temperature too low for molecule formation. 

I Furthermore, although low-velocity shocks are less efficient 

than moderate velocity shocks, they are more numerous 

and depending on the statistical distribution of shocks they 
could potentially account for a higher fraction of the exci- 
tation and formation of molecules. 

The presence of shocks in the line of sight also has im- 
portant consequences on the emission of C + . Indeed, C + is 
usually considered as a good tracer of the ambient UV ra- 
diation field, and indeed the mere presence of the C + ion is 
conditioned to its photo-ionisation. But we show here that 
the emission of C + in both the SQ and Chamaeleon has a 
significant contribution from gas heated by the dissipation 
of kinetic energy rather than by UV photons as is the case 
in photon-dominated regions (PDR). 

The SQ H2 excitation diagram was fit with low and 
moderate velocity C-shocks in dense UV-shielded gas (nH = 
10 4 cm" 3 , G = 0) in earlier studies (jGuillard et al. 2009). 
We show that a better fit can be obtained with diffuse irra- 
diated gas (nn = 10 2 cm~ 3 , Go = 1) including J-shocks. If 
this second solution is the appropriate one, C + is a signif- 
icant coolant of the shocked medium. This prediction can 
be tested with the Herschel observatory. The presence of 
the high-velocity J-shocks can also be tested through the 
emission in HVlines coming from rovibrational levels. 

01, H2O, OH and CO provide additional diagnostic lines 
which could help constraining the physical state of the gas 
and the dissipation processes. OI is seen to be significant at 
larger densities whereas C + is significant for lower densities. 
Molecules such as H2O, OH and CO evacuate less than a 
percent of the total kinetic power except at large densities 
where H2O in C-shocks can contribute to as much as 10%. 

The remainder of the incoming power in our shocks is 
converted into a flux of magnetic energy. The differences in 
total pressure (magnetic plus thermal) between the shocked 
regions and the regions which have not been shocked will 
initiate new motions: the energy stored in the post-shock 
compression will eventually give rise to expansion motions, 
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with an effective transfer back to kinetic energy. Thus, 
shocks are actors which contribute to the equipartition be- 
tween various forms of energy, thanks to a continuous re- 
cycling of the gas throughout successive shock structures. 

The main caveat of our study is that we assume the 
relaxation of the post-shock pressure takes place instanta- 
neously: namely, our model accounts for the pre-shock gas 
and the shocked gas but does not incorporate the expansion 
of the gas after the shock. We shall attempt to model this 
phase in future studies. Another caveat of our study 
lies in the implicit assumption that all shocks are 
steady. Indeed some shock velocities are prone to 
instabilities and even quasi-steady models are in- 
valid in this case (see [Lesaffre et al.l l2004bl for ex- 
ample) . Less important caveats are the shock orientation, 
its curvature, the neglect of intermediate ages of shocks (for 
ages lower than about 10 4 yr, we should consider CJ-type 
shocks or non-steady shocks). We believe that these minor 
caveats are washed out by the fact that we fit a statistical 
collection of shock velocities. Final ly we do not treat UV 
pumping in the H2-population but iMonteiro et al.l (|l988h 
showed this starts to make a difference only for H2 levels 
above J = 5 and for Go 3> 1. 

Our PDF fitting technique might prove useful to ac- 
count for single shocks with more complicated geometries. 
For example, a single bow-shock effectively encompasses a 
continuous collection of velocities depending on the angle 
at which each fluid parcel of pre-shock gas impinges on the 
bow-shock. In such a framework, the shape of the lines is 
directly related to the dynamics inside the shock, and not 
simply due to the relative motions between each shocks as 
in the statistical distribution of shocks we infer in the SQ. 
In such cases, it will therefore be interesting to predict and 
use the line shapes as additional observational constraints 
on the shock physics. 

Finally, as we already mentioned, the dissipation of tur- 
bulent motions does not occur only in shocks. It will be 
interesting in future work to assess on the one hand what 
fraction of the energy dissipation takes place in vortices, 
current sheets or shocks, for example. On the other hand, it 
will b e worth comparing vortex models as in lGodard et al.l 
(2009) to our own shock models in order to check whether 
we can disentangle observationally the respective signatures 
of shocks and vortices. 
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